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■ This work is the continuation of the earlier efforts to apply the 

mean field approximation to the continuum world sheet formulation 
of planar </> 3 theory. The previous attempts were either simple but 
without solid foundation or better founded but excessively compli- 
cated. In this article, we present an approach both simple, and also 
systematic and well founded. We are able to carry through the leading 
order mean field calculation analytically, and with a suitable tuning 
of the coupling constant, we find string formation. 
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1. Introduction 

Over the last couple of years, the present author and Charles Thorn 
have been pursuing a program of summation of planar graphs in field the- 
ory [1,2,3,4]. Because of its simplicity, the field theory that has been most 
intensively investigated so far is the 3 theory, although Thorn and collab- 
orators have made considerable progress in extending the program to more 
physical theories [5,6,7]. The basic idea, which goes back to 't Hooft [8], is 
to represent Feynman graphs on the world sheet by a suitable choice of the 
light cone coordinates, 't Hooft 's representation, which was non-local, was 
later reformulated as a local field theory on the world sheet by introducing 
additional non-dynamical fields [1]. This reformulation has no new physical 
content; it merely reproduces perturbation theory. However, it provides a 
setup well suited for the study of string formation in field theory. This is an 
old problem that has attracted recent renewed interest [9,10,11,12] following 
the discovery of AdS/CFT correspondence [13,14]. 

Our approach to the problem of string formation starts with the world 
sheet description of the 3 field theory mentioned above, and we look for 
the phenomenon of "condensation" of Feynman graphs. This phenomenon 
will be defined more precisely later in the paper, but roughly it means that 
the lines that form Feynman graphs on the world sheet become dense, and 
graphs of asymptotically infinite order dominate the perturbative expansion. 
Furthermore, the original non-dynamical world sheet fields become dynami- 
cal and string formation takes place. Whether the scenario described above 
really happens is of course a question of dynamics. So far, the only tool used 
to investigate this problem in the present context has been the mean field 
or the self consistent field approximation [2,3,4]. The accuracy of the mean 
field approximation is questionable; however, one can hope that at least a 
qualitative understanding of the relevant dynamical issues would emerge. 

The main virtue of the mean field method is its simplicity. There is, 
however, no unique way to do the mean field calculation, and it all depends 
on the choice of the order parameter. In reference [2] , the simplest choice was 
made for the order parameter by taking it to be the expectation value of a 
local field. This makes the subsequent calculation quite tractable. However, 
this early attempt, at least in its continuum version, relied on a number 
of questionable steps and approximations that are hard to justify. In order 
to overcome these difficulties, in references [3] and [4] the order parameter 
was taken to be the expectation value of two fields at different points (two 
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point function). This also has the advantage of providing a better probe 
of the problem; but the disadvantage is that the calculation becomes too 
complicated to carry out analytically. 

In this article, we revisit the earlier calculation [2] based on a simple order 
parameter to see whether the difficulties associated with it can be overcome 
without sacrificing its basic simplicity. We will mainly focus on the treatment 
based on the continuum world sheet given in this reference. Reexamining this 
treatment, we identify the following problems: 

a) The boundary conditions on the world sheet were imposed only approxi- 
mately through the so called (3 trick. The exact boundary conditions could 
only be recovered in the problematic (5 — > oo limit. 

b) If we try to impose the boundary conditions exactly by means of Lagrange 
multipliers, as was done in [1] and as we shall do here, we avoid the problem 
discussed in a), but instead we encounter another difficulty The action has 
then an unfixed gauge invariance, which can lead to ill defined results. 

c) There are two kinds of fields in the problem: The matter and the ghost 
fields. The contribution of each sector is quadratically divergent, but there 
is a subtle partial cancellation between them. Unless great care is exercised, 
the result can depend on the regulation scheme used. 

d) The use of light cone coordinates obscures the covariance of the theory. An 
approximation scheme, such as the mean field approximation, could easily 
violate Lorentz invariance. 

In this article, we propose the following improvements to overcome the 
problems listed above: 

a) The boundary conditions will be imposed exactly by means of Lagrange 
multipliers. 

b) To avoid the resulting gauge invariance problem, we introduce a gauge 
fixing term. This is probably the most important new idea introduced in the 
present work. As is well known, working with a gauge unfixed action can 
lead to lots of confusion. 

c) To keep track of the cancellation between the mattter and ghost fields, 
we impose supersymmetry on the world sheet. This idea was already in 
the air in [1], what we have done here is to formulate it explicitly. Also, 
in treating divergent quantities, we first combine the contributions of the 
matter and ghost sectors, and then we regulate the sum in a way that does 
not spoil the cancellation between them. The answer obtained in this fashion 
is unambiguous. 

d) There is a particular subgroup of the Lorentz group under which the 
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light cone variables transform linearly. In particular, under boosts along the 
1 direction, the variables x ± and p^ 1 scale (see section 2 for the definition 
of these variables). All of this is familiar from the light cone treatment of 
the bosonic string, where, among other things, the importance of invariance 
under this special boost was recognized [15]. In the present context, this was 
discussed in reference [4]. In this paper, following the treatment given in [4], 
we will try to preserve invariance under this special subgroup at each step. 
We will not, however, try to investigate invariance under the full Lorentz 
group. In string theory, full Lorentz invariance in the light cone framework 
is realized only for a fixed critical value of D, where D is the dimension of 
the transverse space [15]. If the same is true here, this problem is clearly 
beyond the scope of a leading large D approximation. 

When these improvements are incorporated in the mean field approxima- 
tion of reference [2], the result is a systematic approximation scheme without 
any of the ambiguities encountered in [2]. It is also simple enough so that 
we are able to carry out the leading order calculation analytically. With a 
suitable tuning of the coupling constant, we find string formation, and in ad- 
dition, we discover that a new string mode has been dynamically generated. 

In organizing this paper, the goal was to present a complete and self- 
contained treatment which should be intelligible even to a reader unfamiliar 
with the previous work on the subject. As a result, there a good deal of 
overlap with reference [2], and some overlap with references [3] and [4]. When 
we preview the rest of the paper below, we will try to make clear what is 
new and what is a review or eloboration of the material covered in the earlier 
work cited above. 

In section 2, we briefly review both the rules for Feynman graphs in light 
cone variables [8] and the local field theory on the world sheet from which 
these graphs follow [1]. We also discuss the transformation properties of the 
fields under the special boost mentioned earlier, which manifests itself as a 
scale transformation on the world sheet. This is an abridged version of a 
more complete discussion given in [4]. 

In section 3, world sheet supersymmetry is introduced and a supersym- 
metric free action Sq is constructed. This is a new idea in the present context. 
5*0 differs from the corresponding action introduced in [1] and used in all the 
subsequent work in the structure of the terms involving ghosts, and also in 
the presence of a relatively insignificant bosonic field r required by SUSY. 
The cancellation of quantum corrections between matter and ghost fields, 
which was the reason for the introduction of ghosts in the first place, is now 
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guaranteed by supersymmetry. 

The boundary conditions accompanying So are enforced by a term in the 
action, Si, given in section 4. Both the boundary conditions and Si are 
substantially the same as those in [1]. To express the new term in the action 
fully in the language of field theory, one needs world sheet fermions which 
were introduced earlier [1,4]. We have found it slightly more convenient 
to work with a somewhat different set of fermions, although it is easy to 
show by means of a Bogoliubov transformation that the two are completely 
equivalent. 

In section 5, we show that the action constructed so far is invariant under 
a simple gauge transformation. It is therefore important to fix this gauge, 
and we show how to do it. We also note that it is possible to introduce 
some arbitrary constants in the boundary and gauge fixing terms in the 
action. Although the exact theory is not sensitive to these constants, we 
keep them around to see how the affect an approximate calculation. The 
material covered in this section is completely new. 

The mean field approximation which is at the basis of the present work is 
discussed in section 6 from the point of view of the large D limit. All of this is 
standard material, well known from the solution of the vector models in in the 
large iV (in this case, large D) limit [16]. The only thing new here compared 
to [2] and [3] is the manner in which the singular determinants resulting 
from integration over the matter and ghost fields are regulated: The two 
determinants are combined into an expression less singular than each one 
seperately, and regulating the combined expression, we get an unambiguous 
and scale invariant answer. 

In section 7, an effective action is constructed and evaluated by the saddle 
point method in the large D limit. This effective action is pretty close to but 
still different in detail from the one derived in [2] . A standard bosonic string 
action with positive slope emerges from this calculation. The important 
question is then whether this slope is finite. We find that, by suitably tuning 
the coupling constant of the 3 theory, the slope can be made finite. This 
tuning can be regarded as renormalization: A cutoff dependent coupling 
constant is traded for the finite slope parameter. 

It seems somewhat surprising that starting with an unstable field theory, 
no sign of instability has so far appeared in the string picture. One possibility 
is that we have not gone far enough. The calculation of the intercept, which 
we have not undertaken here, may show that, as in the bosonic string, some 
lowest lying states are tachyonic. Or, it may be that the instability is not 
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visible in the leading order of the large D limit. These possibilities are 
discussed in at the end of section 7. 

So far, all the calculations were carried out in the leading order of the 
large D limit. In section 8, we compute a non-leading correction in this 
limit by expanding the composite field p (see section 4 for its definition) 
around its mean value p . In addition, p is assumed to be slowly varying, 
and an expansion up to second order in powers of derivatives of this field is 
carried out. This is essentially a repetition of the computation done in the 
Appendix B of [2] from the standpoint of the present approach. We find that, 
from the world sheet perspective, p becomes a dynamical, propagating field, 
and from the string perspective, the string acquires an additional mode, with 
the same slope as all the other modes. Finally, in the Appendix, we show 
that, the mean field computation presented here is completely equivalent to 
the standard large N (D) treatment of vector models. 

2. A Brief Review 

The Feynman graphs of massles 3 have a particularly simple form in the 
mixed light cone representation of 't Hooft [8]. In this representation, the 
evolution parameter is x + also denoted by r, and the conjugate Hamiltonian 
is p~, and the Minkowski evolution operator is given by 

T = exp{-ix + p-). (1) 

In this notation, a Minkowski vector v M has the light cone components 
(v + ,v~,v), where v ± = (v° ± v 1 )/^, and the boldface letters label the 
transverse directions. A propagator that carries momentum p is pictured as 
a horizontal strip of width p + and length r = x + on the world sheet, bounded 
by two solid lines (Fig.l). The lines forming the boundary carry transverse 
momenta qi and q 2 , where 

P = qi - q 2 , 

and the corresponding propagator is given by 

More complicated graphs consist of several horizontal solid line segments 
(Fig. 2). The beginning and the end of each segment is where the 3 interac- 
tion takes place, and a factor of g is associated with each such point, where 
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Figure 1: The Propagator 



g is the coupling constant. Finally, one has to integrate over the position of 
the interaction vertices, as well as the momenta carried by the solid lines. 
We note that momentum conservation is automatic in this formulation. A 
typical light cone graph is pictured in Fig. 2. 

It was shown in [1] that the light cone Feynman rules described above can 
be reproduced by a local field theory on the world sheet. The world sheet 
is parametrized by coordinates o along the p + direction and r along the x + 
direction, and the transverse momentum q is promoted to a bosonic field 
q(<7, r) defined everywhere on the world sheet. In addition, two fermionic 
fields (ghosts) b(a, r) and c(a, r) are needed. In contrast to q, which has D 
components, b and c each have D/2 components, where D is the dimension 
of the transverse space, assumed to be even. The action on the world sheet, 
with the Minkowski signature (+,-), is given by 

So = £ + da jT' dr (V • d - iq' 2 ) , (3) 
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Figure 2: A Typical Graph 

where the prime denotes derivative with respect to a. This action is supple- 
mented by the Dirichlet conditions 

q = 0, b = c = 0, (4) 

on the solid lines (boundaries), where the dot denotes derivative with respect 
to t. Imposing Dirichlet conditions on q is equivalent to fixing them to be r 
independent on the solid lines and then integrating over them. If we fix q to 
be qi and q2 on two adjacent solid lines and solve the equations of motion 

q" = 

subject to these boundary conditions, we find that the the resulting classical 
action reproduces the exponential factor in eq.(2). However, there is also an 
unwanted quantum contribution, given by 

-\Ddet(dl) 

which is exactly cancelled by the corresponding determinant resulting from 
integrating over b and c. 



7 



The action formulation described above, which was extensively used in 
the previous work [2,3,4]. Although it is basically correct, it has some unsat- 
isfactory features. For one thing, although we have nothing new to say about 
this problem in this paper, the factor of 1/ (2p + ) in front of the exponential in 
eq.(2) is missing. Also, the splitting of the ghost fields into two components 
b and c is somewhat unnatural; it leads to the artificial condition that D is 
even, and it does not look rotationally invariant, although there is of course 
no real violation of rotation symmetry in the transverse space. In the next 
section, we will show that the introduction of supersymmetry on the world 
sheet leads to a more symmetric ghost sector and the condition that D be 
even is no longer needed. We also think that it results in a more natural and 
elegant approach. 

Finally, we would like to discuss briefly the question of Lorentz invariance. 
This is a non-trivial problem, since the use of the light cone coordinates 
obscures the Lorentz transformation properties of the dynamical variables. 
There is, however, a special subgroup of the Lorentz group under which 
the light cone coordinates have simple linear transformation properties 3 . If 
L i j are the angular momenta and Ki are the boosts, the generators of this 
subgroup are 

Li j, M+- = K U M +>i = K i + L lii , (5) 

where indices % and j run from 2 to D + 2. It turns out that invariance under 
all the generators except for one is rather trivially satisfied by the propagator 
(2) or the field theories (3). The non-trivial transformation, generated by Ki, 
corresponds to scaling of x + and p + by a constant u: 

x + — > x + /u, p + ^p + /u, (6) 

and the tranverse momenta q are unchanged. The reason this transformation 
is critical is that although it is easy to construct classically scale invariant 
theories, this symmetry is in general broken by quantum corrections. This is 
familiar from the study of conformal invariance in field theory. 

Now, let us examine the scale invariance of the action (3) [4]. Since the 
coordinates a and r must transform like p + and x + , the transformation laws 
of the fields are 

q(<7, r) — > q(ii<7, ut), b,c(a,r) — > b,c(ua,ur), (7) 

3 In this article, we will not investigate invariance under the remaining Lorentz gener- 
ators. If the string analogy is to be trusted, this is where the critical dimension becomes 
important [15], and a large D approximation is clearly inadequate for treating this problem 
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and the classical action is invariant under this transformation. On the other 
hand, the quantum corrections are singular and need a cutoff for their defini- 
tion. This cutoff would break scale invariance, were it not for the cancellation 
between the ghost and matter fields. So what we have here is a potential 
violation of scale invariance by an anomaly, which is eliminated by the can- 
cellation between the matter and ghost sectors. This is nothing but the 
cancellation of the determinants discussed earlier. We would like to stress 
that the quantum contribution in question is a world sheet effect; it poten- 
tially present even in the case of a free propagator and it has nothing to do 
with the target space ultraviolet divergences. 

In addition to the scaling of the bulk, we have to consider the scaling be- 
haviour of the boundary conditions given by eq.(4). These are scale invariant 
as they stand, but the integration over the position of the boundaries is not 
invariant, since the a coordinate scales. The factor of l/(2p + ) provides the 
measure needed to make this integration scale invariant. Although we will 
not present here a general recipe for the inclusion of this factor, we will be 
careful to preserve the scale invariance of various integrals that occur in the 
course of the mean field calculation. It can be shown that, in any case, this 
factor does not contribute to the leading order of the mean field calculation 



where 6^2 are the usual anticommuting coordinates and all the fields repre- 
sented by boldface letters are vectors in the D dimensional transverse space. 
The two SUSY generators are(i = 1, 2) 



[2,3]. 



3. SUSY On The World Sheet 



We generalize the momentum q to form a SUSY multiplet: 



Q = q + 1 b + 2 c + 1 2 r, 



(8) 




(9) 




= 0. 



(10) 
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where ei j2 are anticommuting infinitesimal parameters. We also define two 
covariant derivatives 

A = J- - Oid a , (12) 

which anticommute with the SUSY generators and satisfy the same equations 
as (10), apart from the change of the sign of d c . 

Having introduced the framework of supersymmetry, we write down the 
supersymmetric analogue of the action (3): 

S -> £ da J dr J d 2 9 (~^iQ • AQ) 

= £ + da J dr ("^(q' 2 + b • b' + c • c' + r 2 )) , (13) 

supplemented by the boundary conditions 

q = 0, b = 0, c = 0, r = 0, (14) 

on the solid lines. We now have twice as many ghost fields as before, but they 
appear with only first order derivative in a, and as a result, there is again 
complete cancellation between the determinants( quantum corrections) com- 
ing from the integration of the matter and ghost fields. It is clear that this 
cancellation is a consequence of supersymmetry. In fact, the main motivation 
for introducing SUSY was to obtain this cancellation as a systematic conse- 
qence of a symmetry, and not as some accident. In the rest of this article, 
we will exclusively use this supersymmetric form of the action. 

Although there is supersymmetry in the bulk of the worldsheet, it is bro- 
ken by the boundary conditions, since the condition on q differs from those on 
b and c. This breakdown of SUSY is essential, since supersymmetric bound- 
ary conditions would lead to a complete cancellation between the matter and 
the ghost sectors, resulting in a trivial propagator. This has no effect on the 
cancellation of the quantum contributions, since the cancellation occurs in 
the bulk and it is insensitive to the boundary conditions. 

Since we now have a new expression for So, we have to reinvestigate 
the invariance under the scaling transformations (6,7). The action given by 
eq.(13) is invariant if the fields transform as 

q(<7, r) — > <\(ua, ut), b(<7, r) — > \fuh(ua, ut), 

c(a,r) — > y/uc(ua, ut), r(a, r) — > uv{ua, ut). (15) 
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Again, a potential quantum anomaly that could violate scale invariance is 
cancelled as a consequence of supersymmetry. 

4. The World Sheet Action 

In this section, we review the construction of an action that incorporates 
the boundary conditions (14), which will be implemented by introducing a 
bosonic Lagrange multipliers y(er, r) and z(er, r), and the fermionic Lagrange 
multipliers b(<j, r) and c(a, r) [1]. The corresponding term in the action is 

Si = J* da J d,Tp(y-q + b-b + c-c + z- t) , (16) 

where p(a, r) is equal to 1 on the solid lines and it is equal to elsewhere. 
This definition is singular on the continuum world sheet; we will give a more 
precise definition on the discretized world sheet below. The factor p ensures 
that the boundary conditions are imposed only on the solid lines; elsewhere, 
p vanishes and there are no constraints on the fields. 

In the earlier work [1,2,4], it was shown how to construct p in terms of 
a fermionic field. We will present here a slightly different version of this 
construction. In order to avoid singular expressions, it is best to start by 
discretizing the a coordinate into segments of length Aa = a, with p + = Na. 
The parameter a plays the role of a cutoff, which will in any case be needed 
later on. The specific form of this cutoff is not important; for example, a 
cutoff in the Fourier modes conjugate to a would serve just as well. Fig. 3 
shows N equally spaced lines parallel to the r direction. For the convenience 
of exposition, the absence of a solid line in a given position is pictured by 
the presence of a dotted line in the same position. In other words, the solid 
lines mark the boundaries and the dotted lines fill the bulk of the world 
sheet. On each line, we introduce two fermionic variables ipi{a ni r) and their 
conjugate ipi(cr n , r), where i — 1, 2, a n — na, and n is an integer in the range 
< n < N . The free fermionic action is given by 

S 2 = Yl / dri^((T n ,T)ij(a n ,T), (17) 

n J 

and the fermions satisfy the usual anticommutation relations: 

[$i(a m , r),i}j(a n , r)]+ = 5 i:j 5 m>n . 
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Figure 3: Solid And Dotted Lines 



The function of the fermions is to keep track of the solid and dotted lines. 
The vacuum state, defined by 

^|0> = 0, 

corresponds to the trivial situation where all the lines in the graph are dotted. 
The r independent state corresponding to a single solid line at a = a n is 

\<?n) = 4>l(cr n ,T)4> 2 ((Tn,T)\0). (18) 

This is an eternal solid line, extending indefinitely for both positive and 
negative r, since the operators ip, and hence the above state, are independent 
of t by the equation of motion following from (17). Several solid lines are 
represented by the state 

n^K)^K)|o). 

n 

We note that, a state with a double solid line, which has no graphical inter- 
pretation, vanishes by fermi statistics. Having set up our fermionic system, 
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we can express p in eq.(16) in terms of fermions: 




) = ;E ipi(a n ,T)ipi(a n ,T), 



(19) 



i=l,2 



and it is easy to check that p is 1 if there is a solid line located at a = a n , 
and it is zero if the line is dotted. The ground state expectation value of 
this composite field, po, will play an important role in what follows. In 
any finite order of perturbation theory, where the density of solid lines is 
essentially zero, po is zero. On the other hand, a non-zero po means that 
a portion of the world sheet with finite area is occupied by the solid lines, 
which we interpret as condensation of solid lines. It is also clear that this 
can only happen if large (infinite) order graphs dominate the perturbation 
series. It will be shown later that, at least in the mean field approximation, 
the condensation of solid lines leads to string formation. One can think of po 
as an order parameter that distinguishes between the stringy phase and the 
perturbative phase of the same theory. 

So far, all the lines, whether solid or dotted, are eternal. We need an 
interaction term in the action which will convert dotted lines into solid lines 
and vice versa. Remembering that the transition between solid and dotted 
lines is accompanied by a factor of the coupling constant g, we set 



and it is easy to show that this term in the action does the required job. 

As we mentioned earlier, the fermions introduced in this paper are some- 
what different from those used in the earlier work [1,2,4]. However, it is not 
difficult to show that the two are connected by the Bogoliubov transformation 



and, as a result, they are physically equivalent. 

Finally, we will consider the continuum limit, with a — > 0. The dictionary 
for taking this limit is 



» 



(20) 



n 



V>1 <^> Ipl, 



E 



n 




1 




1 



(21) 



a 
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For the continuum fermions, we use the same notation as the discrete ones, 
but with a n replaced by a. They satisfy anticommutation relations similar to 
the discretized version, with the Kroenecker delta replaced by the Dirac delta 
function. The expressions for S2 and S3 (eqs.(17) and (20)) remain unchanged 
if one replaces the sum over a n by integration over a. In particular, there 
is no explicit dependence on a. The state corresponding to a solid line at 
a = a' is now represented by 

W) = ^1(^,^)^2(^,^)10), 

and 

p(a,r)\a')=5(a-a')\a'). (22) 

In what follows, although we will mostly work with the continuum fermions, 
from time to time we will also need the world sheet with discretized a to 
have well defined expressions. 

The total action is the sum of various pieces given by eqs.(13),(17) and 
(20): 

S = S + S 1 + S 2 + S3. (23) 

As it stands, this action suffers from a serious problem: It is not well defined. 
We will discuss this problem and present its solution in the next section. 

We close with a discussion of the scaling properties of the fermions. In- 
variance of the canonical algebra, or of the free action for the continuum 
fermions requires the transformation law 

i/j(<j,t) — > y/uip(ua,ur), ^(<j, r) — > \/u?P(u(J,ut), (24) 

and as result 

p(<7, r) — > u p(uct,ut), P~(<J,t) — > u p(ua,UT). (25) 

We also note that, from its definition Na = p + , it follows that the cutoff 
parameter a must scale as p + : 

a — > a/u. (26) 

It is now easy to check that, with the possible exception of the interaction 
term S3, all the terms of the fermionic action are scaling invariant. The 
interaction term violates scaling, unless the coupling constant g is allowed 
to transform. Of course, the original field theory coupling constant is a 
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Lorentz scalar and cannot transform. However, the coupling constant on the 
world sheet is closely related to but not identical to the field theory coupling 
constant, and it need not be a Lorentz scalar. In fact, Lorentz invariance 
requires that 

g -> ug 

under scaling. The simplest way to secure this is to use p + , the only other 
physical parameter at our disposal and make the replacement 

« - < 27 > 

p + takes care of the scaling (see eq.(6)), and the newly defined g is a true 
scalar. 

5. The Gauge Fixed Worldsheet Action 

Consider eq.(16) for Si. Since p vanishes everywhere except on the solid 
lines (the boundary), the integrand is independent of the Lagrange multipliers 
y, b, c and z in the bulk of the world sheet. As a result, the functional 
integration over the bosonic variables is divergent 4 and the corresponding 
integration over the fermionic variables vanishes. The result is an ill defined 
expression of the type infinity times zero. This is similar to what happens in 
gauge theories before gauge fixing. In fact, Si is invariant under the gauge 
transformation 

y -»• y + pyo, b^b + pb , 

z — > z + pz , c^c + pco, (28) 
where yo, z , b and c are arbitrary functions and p is defined by 

for discretized a and 

P=--P = \Y.^ ( 29 ) 

for continuous a, where p is given by (19) in the first case and by (21) in the 
second case. In either case, since 

P P = 0, (30) 

4 I thank Charles Thorn for stressing this point, although at the time I did not think it 
was important. 
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invarinace of Si under the transformations (28) follows. It is this gauge 
invariance that is responsible for the singular behaviour of the functional 
integrals over the Lagrange multipliers mentioned above. 

The cure for this problem is to gauge fix the action, but only in the bulk 
of the world sheet (on the dotted lines), where p = 0. We will also demand 
the gauge fixing term to be supersymmetric, since our philosophy is to keep 
intact the SUSY in the bulk and violate it only on the boundaries. We 
therefore promote the Lagrange multipliers into a supersymmetric multiplet 

Y = z + 0ic + 2 b + 0i0 2 y, (31) 
and write the gauge fixing part of the action as 

S A = £ da J dr J d 2 9^apY Q 

= J da J dr ap (y • z + b • c) , (32) 

where a is a gauge fixing parameter. This is the simplest gauge fixing term 
which is supersymmetric and which vanishes on the boundaries (solid lines) 
where p = and no gauge fixing is needed. The previously singular functional 
integral over the bulk of the world sheet where p = is now equal to unity: 



Vy VzVbVc exp(iS 4 ) = 1. (33) 

We note that the a dependent integration measures of the bosonic and 
fermionic functional integrals cancel. We believe that the absence of gauge 
fixing on the continuum world sheet caused some of the problems encountered 
in the earlier work. 

There is one more technical issue which we have to address here. There 
is some arbitrariness in the equation for (16) for Si; it could be replaced by 
the following more general expression: 

Si -> jT da j dr p (ft y • q + ft b • b + ft c • c + & z • r) . (34) 

Here, /3i,2,3,4 are arbitrary constants. They can be eliminated by absorbing 
them into the definition of the Lagrange multipliers. So long as one is dealing 
with the exact expression for the action, the introduction of the /3's changes 
nothing, and one could just as well set them all equal to one, as in eq.(16). 



16 



However, if an approximation scheme is used, the results may well depend 
on these constants. These remarks also apply to the gauge fixing parameter 
a: The exact theory is independent of this constant but an approximate 
calculation may introduce some dependence. In the next section, when we 
carry out a mean field calculation, we will be able to see to what extend our 
results are sensitive to the choice of these constants. 

Finally, we have to make sure that the gauge fixing term S4 and also S\ 
are scale invariant. This is indeed the case for S4 if under scaling 

y(o-,r) -> y(ua,ur), b(a, r) -> y/ub(ua, ut), 

c(<7,t) — > ^/uc{ua, ut), z((T,r) ->«z(tt(j,Mr), (35) 

and eq.(25) is taken into consideration. Now that we know the scaling prop- 
erties of all the fields, we can check Si (eq.(16)). The first three terms are 
indeed invariant as they stand, but in the last term, we have to let 

6. The Mean Field Approximation 

In the last section, the final form of the action which is supersymmet- 
ric, gauge fixed and scale invariant was worked out. Since the mean field 
approximation will be based on it, we start by writing it down in full: 

n=4 

s = £s n 

n=0 

= J* da J dr( - i(q /2 + b • b' + c • c' + r 2 ) + p{^ y • q + (3 2 b ■ b 

+ (3 3 c ■ c + [3 A p + z- r) + ap (y • z + b ■ c) + iiptp + ^(^1^2 + ^2^1)) > 

(36) 

where p and p are given by (21) and (29). 

The mean field approximation was developed and applied to the world 
sheet action in the earlier work [2,3,4]. Here, we need the simplest version of 
it used in [2]. We notice that eq.(36) represents a vector model, which can be 
solved in the large D limit [16]. The standard approach is to replace scalar 
products of various vector fields, such as the bilinear terms in the expression 
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for So (eq.(13)), by their vacuum expectation values, which are then treated 
as classical c-number fields. The functional integral over the remaining fields 
is carried out exactly, and the resulting effective action is minimized with re- 
spect to the classical fields. This method is known to be equivalent to a large 
N ( in this case, D replaces N) saddle point calculation [16]. Instead of the 
approach sketched above, we find it much simpler to replace the bilinears in 
the fermions, p and p by their classical expectation values. In the Appendix, 
we show that this is completely equivalent to replacing the bilinears in vector 
fields by their expectation values. 

The problem is considerably simplified by starting with the setup where 
the total transverse momentum p carried by the whole graph is zero: 

p= f P daq=0. (37) 
Jo 

This configuration can always be reached by a suitable Lorentz transforma- 
tion. It allows us to impose the periodic boundary conditions 

q( ( 7 = 0,r)=q(a = p + ,r). (38) 

This setup is translationally invariant in both the a and the r directions, 
and we shall see that translation invariance will play an important role in 
simplifying the mean field calculation. 

We start by explicitly introducing the composite field p by adding a new 
term AS to the action: 

S -> S + AS, 
AS = J da J dm 

where k is a Lagrange multiplier, and p is given in terms of p through eq.(29). 
In the large D limit, we can treat k and p as classical fields (See the Ap- 
pendix). In other words, we make the replacement 

K -y k = (k), p -> po = (p), p -> po = (p), (40) 

where () denotes the expectation value in the ground state of the field in 
question. Translation invariance on the world sheet means that both kq 
and po are constants independent of the coordinates a and p. With this 
simplication, it is possible to carry out explicitly all the functional integration 



2 2 



(39) 
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over the fields. We first consider the integration over the fermions; that part 
of the action involving the fermions is given by 

Sf = Io ^/ rfr (^ + ^(^ + «i) + yE^- (41) 

Instead of working with the action, we find it more convenient to diagonalize 
the corresponding Hamiltonian. In order to avoid singular expressions, we 
first discretize the a coordinate as in eqs.(17) and (20). There is a complete 
decoupling of the dynamics in the a direction; as a result, the total Hamil- 
tonian can be written as a sum of iV mutually commuting Hamiltonians: 

H = Y, H m (42) 



where, 



and 



H n = ~ WiM* + iMi) + y E^j , (43) 



9 = g/p + - 

We observe that H n acts on the two states 

|0> 

corresponding to a dotted line at a = a n and 

\n) = ^i((7„)^ 2 (o"n)|0) 
corresponding to a solid line at the same position as a two by two matrix: 

H n \0) = -g'\n), 

H n \n) = -<7'|0)-/co|n). (44) 

The eigenvalues are 

Ef = ~Ko±y^ + Ag», (45) 

where the square root is defined to be positive. Since we are interested in 
the ground state, we have to pick the lower energy Ej at each a = a n and 
add up to get the total fermionic energy: 

Ef = NEj = V —E r (46) 
19 



Next, we will carry out the functional integrations over the vector fields 
in eq.(36), setting 

1 

P~> Po, P -> Po = Po- 
tt 

We first split the action into S m , the matter action, and S g , the ghost action: 

S m = J o da J dr (-^(q' 2 + r 2 ) + Po(PiY • q + /3 4 p + z • r) + ap y ■ z^j , 

S 9 = £ da J dr (-^(b • b' + c V) + p (/3 2 b + /3 3 c • c) + ap b • . 

(47) 

The integrations over y, r, b and c are Gaussian and they can be done 
trivially. The (singular) Jacobians coming from integrations over the matter 
and the ghost fields cancel, with the result, 



ap 



(48) 



We note that S m is the(Minkowski) world sheet action for a string, with 
the slope parameter a', where 

Aa-=( ^ P+p2 A 2 . (49) 
V "Po / 

It may seem like string formation is almost automatic; however, the string 
picture breaks down if the slope is zero , which happens for p = 0. The 
parameter po is therefore the order parameter that distinguishes between 
the stringy and the perturbative phases of the same field theory. Roughly 
speaking, since p measures the fraction of the world sheet area occupied 
by the solid lines, each graph in perturbation theory corresponds to po = 0. 
This is because any finite number of solid lines, being one dimensional, have 
vanishing area, and as to be expected, perturbative field theory is then the 
zero slope limit of the string theory. A non- zero slope requires po 7^ 0, 
which means that the solid lines condense to occupy a finite fraction of the 
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area of the world sheet. Therefore, po serves as an order parameter which 
distinguishes between two phases: po — in the perturbative field theory 
phase and po 7^ in the stringy phase. In the next section, we will find that 
Po 7^ in the mean field approximation. 

Since S mt9 have quadratic dependence on the fields, the functional inte- 
grals can be carried out. Defining 

exp(iS m>g ) = exp(iS^ g ), 

we have, after Euclidean rotation in r, 

S e m = ~DTr\n(-dl-A^) : 
S e g = DTr\n(-dl-B 2 dl), (50) 

where, 

A = /Wp8 B = . Ml 

Rewriting the Tr ln's in terms of momenta ko and k\ conjugate to r and a 
gives 

S ^Sf„ + ^^/*pn(|pl±^), (5!) 

with 

f,n = 2 ™ 
1 p+ ' 

This is only a formal result, since the integration over ko and summation 
over n lead to a quadratic divergence, and before one can make sense of it, 
it must be regulated. In fact, we are precisely interested in the coefficient 
of this quadratic divergence, which, after it is regulated by a cutoff, will be 
the dominant term in the answer. Since the answer is somewhat sensitive to 
the cutoff procedure, we will try to explain the motivation for the regulator 
we use. First, we make use of a simplification: Only the leading cutoff 
dependent part of the answer is of interest and this dependence is sensitive 
only to the large fc" and large ko behaviour of the integrand in the above 
equation. Therefore, we can safely replace the summation over n by the 
integration over the continuous variable k±. By the same token, we can set 
B = in the integrand: 

S e - (r f -n)^r[ dko I dh In (j^r^) • (52) 



8tt 2 J U J 1 \kl + A 2 kl 
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In eq.(51), the contributions of the matter and ghost sectors were com- 
bined into a single log. This was intentional: We are going to regulate the 
combined contributions of the two sectors, rather then regulate them seper- 
ately. To see why, we observe that 

a) The combined term is less singular then each term treated seperately. In 
fact, for a fixed k , the integral over k\ converges, so that we need to regulate 
only the integral over fc . This is no accident: it can be traced back to the 
cancellation of the singular determinants between the matter and the ghost 
sectors. Regulating each term seperately could spoil this cancellation. 

b) We saw in section 2 that scale invariance on the world sheet is necessary 
for Lorentz invariance in the target space. Therefore, the regulated expres- 
sion for S e should be scale invariant. Under scaling, k and hi transform 
as 



and from eq.(50), it is easy to check that A is scale invariant. It then follows 
that the integrand in eq.(52) is scale invariant. This is of course related to the 
cancellation discusssed above, if we recall from section 3 that the cancellation 
of determinants and scale invariance are intimately connected. Again, it was 
essential to combine the matter and the ghost terms to arrive at a scale 
invariant integrand. 

The scale invariance of the integrand in eq.(52) is necessary but not suf- 
ficient for the scale invariance of S e ; one also needs a regulator that respects 
scaling. We have seen that only the ko integration has to be regulated, which 
we regulate by introducing a second cutoff (in addition to a) in the r direc- 
tion. Again the precise form of the regulator is not important, so long as it 
respects scaling. As a simple example, we will consider a sharp cutoff in k : 



where A is the cutoff parameter. The factor of p + is inserted so that the 
limits of ko integration are invariant under scaling. By a change of variables, 
we can evaluate this integral as 



k$ — > u ko, k\ — > u k\ 




(53) 



S e = 



D Tf-Tj 

8tt 2 p+ 



A A, 



(54) 



where, 




(55) 
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We note that 

a) A is a positive cutoff dependent constant. Since this is the only cutoff 
dependent parameter in the result, we may as well replace the original cutoff 
parameter A by A. 

b) The simple linear dependence of S e on A is going to be important in the 
following development. This dependence is a fairly robust result: It follows 
from the change of variable 

k\ — > A hi, 

independent of the details of how the integral over ko is regulated. 

c) The dependence on p + is required by scale invariance, again independent 
of the form of the regulator. 

The preceding discussion leads to the conclusion that S e has the unique 
form given by eq.(54), provided that we combine the matter and ghost de- 
terminants before regularizing and we use a regulator that respects scale 
invariance. 

Finally, we would like to comment on the factor of i in the definition 
of B in eq.(54). If the product P2P3 is real, B 2 will be negative, and from 
eq.(), S e will be complex. This is, of course, a signal for instability. This 
does not concern us here: Since we are interested only in the leading cutoff 
dependence, we use eq.(51), where the B 2 term is absent. However, even for 
the non-leading part of S e , it is possible to avoid this problem by taking the 
product P2P3 to be pure imaginary. We recall that the constants (3 2 and /? 3 
appeared in front of the Lagrange multipliers that set the fields b and c equal 
to zero. If these fields were bosonic, complex values for these constants would 
not be permissible, but for fermionic fields, complex coefficients are allowed. 
Nevertheless, this may still show up as an instability in some non-leading 
order in the large D limit. 

7. String Formation 

We can now put together various terms derived in the last section and 
write down the full effective action S e ^: 



S eff = S e -(r f -T t )(E f +p + K p ) = ^- T ^S, 
87r 2 ap 



S = -^^A + ^(. + ^)- W W 



(56) 
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Now let us go back to the question posed in Section 5: How does the effective 
action depend on the arbitrary constants /3i,2,3,4 that appear in eq.(36) and 
the gauge fixing parameter a? The answer is that, at least within the present 
approximation, it does not depend on /3 2 and (3% at all. Also, the dependence 
on a and j3\ and j3± is rather trivial. Since the cutoff parameter A is arbitrary 
to begin with, by redefining it through 

a 

one can eliminate all the reference to these parameters. 

Next we will search for the saddle point of the effective action in the vari- 
ables po and /to, or equivalently in the variables x and y, and see whether this 
corresponds to the minimum value of the ground state energy. To start with, 
the expression for S can be simplified considerably by a series of redefinitions: 

l-x DA . g DA 

Po = , ^0 = 2— y, 9=—r = —9- (57) 

a p + p + p + 

It is also convenient to define a scale invariant cutoff parameter a' by 

a' = a/p + . 

In terms of these new variables, S is given by 

S = D±F(x,y), (58) 

CI 

where 

2 

F(x, y) = -2xy + y + ^y 2 + g 2 , (59) 

X Ob 

and y ranges from —00 to +00, whereas x is limited to the interval < x < 1. 
A number of features of this expression are worth noting: 

a) There is a factor of D multiplying the whole expression. Therefore, in the 
large D limit, the dominant contribution comes from the saddle point. 

b) Every variable in this expression is scale invariant. 

c) The factor of D appearing in the definition of g is the standard "large N" 
factor [8] needed to have the correct limit. 
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d) Although we have so far introduced two independent cutoff parameters A 
and a, or equivalently, A and a', only the combination A/a' appears in the 
expression for S. 

The integral to be evaluated in the large D limit is 

/+°° /"I 
dy / dx exp(iS eff ), 
-oo JO 

S eff = lF(x,y), (60) 
and the constant 7 is defined by 

A Tf - Ti 

1 = D- . 

a' p + 

First, let us consider the integration over y for a fixed value of x in the 
interval (0,1). This integral can be evaluated exactly in terms of a Bessel 
function, but since we need only the large D result, we will instead use the 
saddle point approximation. The saddle point y s is at 

2x — 1 

d y F(x, y)=0^y s = —j== 9. (61) 

In this equation, both g and the square root are defined to be positive. The 
contour of integration in the complex y plane can be distorted into the curve 
whose equation is 

Re (-2xy + y + ^y 2 + g 2 ^ = 2gVx - x 2 . (62) 

This curve passes through the saddle point, and with the above choice of the 
branch of the square root, as 

Im(y) — > ±00, 

on the curve, 

Im (F(x,y)) — > +00, 

and therefore the integral 



L 



dy exp (i~fF(x,y)) 
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Figure 4: The Integration Contour in the Complex y Plane for g = 1 

is exponentially convergent. Using this contour of integration justifies the 
evaluation of the integral over y by setting y = y s in the integrand. In Fig.4, 
the curve defined by eq.(62) in the complex y plane is pictured for g = 1 and 
x = 1/2, when the saddle point is at y s = 0. The branch cuts of the square 
root run from ig to +ioo and from — ig to — ioo and the contour asymptotes 
the vertical lines Re(y) = ±g. 

After the saddle point evaluation of the integral over y, we are left with 
the integral over x: 

Z — > / dxexp (— -ry f(x)) , (63) 
Jo 
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Figure 5: The Function f(x) for g = 20 



where 



/(*) 



x 



x 



— 2g\/x — x 2 . 



(64) 



The function f(x) is pictured in Fig.5 for g = 20. It has a single minimum 
in the interval 0<a;<lata; = x m ,, which satisfies 



/'W 



9- 



1 - 2X T 



0. 



(65) 



At the minimum, /(x) is negative, corresponding to a negative minimum 
energy 

E = D-^-f(x m ), 
a' p + 

which is the energy of the ground state in this approximation. The situation 
is similar for other positive values of g: There is only a single minimum in 
the interval (0, 1), corresponding to a negative ground state energy. 

We have just seen that the ground state corresponds to a value of x 
between and 1. Let us remember that 



x = — > po = 
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corresponds to the trivial case of a world sheet with all dotted lines. The 
opposite limit of 

x = 1 — > po = 

corresponds to a world sheet with only solid lines. A value of x in between 
these two extremes implies an intermediate world sheet texture: The solid 
and dotted lines each occupy a finite fraction of the area of the world sheet. 
Recalling our earlier discussion following eq.(49), we see that indeed a con- 
densate of the solid lines has formed. 

Let us now see whether a sensible string picture emerges. In particular, 
the slope parameter (eq.(49)) 

a' = ^ P+Pl = f 1 ^ , (66) 
2apo 2a(l — x m )a' 

is the only physical parameter to emerge from the mean field calculation. 
The theory should be renormalized by requiring it to be a finite number 
independent of any cutoff. To see how this happens, we first get rid of the 
irrelevant constants (3^ and a by suitably redefining the cutoff parameter 
a': 

a' -> — Xm . . (67) 

In order to have a finite slope in the limit a' — > 0, x m , and therefore, g should 
also go to zero in the same limit. Solving (65) in the small g limit, we have 

x m » (JJ , (68) 

<*' - i^r (69) 

Therefore, in the limit of the cutoff a' tending to zero, the coupling constant 
g should be fine tuned so that the ratio 

a' 

stays constant. The theory is renormalized by trading the cutoff dependent 
coupling constant g for the cutoff independent physical parameter a' through 
eq.(69). 



and 
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There is one more comparison one can make between the parameters 
of the 3 theory and the string theory: One could try to relate the string 
intercept to the field theory mass and the coupling constant. If we identify the 
mass square of the lowest lying state on the string trajectory with the ground 
state energy corresponding to S (eq.58), the meanfield approximation gives 
a negative cutoff dependent answer. In this article, we will not consider the 
question of renormalization of the ground state energy, which is the same as 
the renormalization of the intercept. There is, however, the following simple 
possibility. Instead of starting with a 3 theory with zero bare mass, we could 
have started with a non-zero cutoff dependent bare mass. Fine tuning this 
bare mass, it may be possible to end up with a finite renormalized ground 
state energy. To carry out this program, however, our formalism has to be 
extended to include a non-zero bare mass. This we leave to future research. 

It may seem surprising that, starting with an unstable field theory, so far 
we have not encountered any sign of instability on the string side. Of course, 
as mentioned above, we have not calculated the renormalized intercept. In 
the end, upon calculating this intercept, just as in the case of the bosonic 
string, some of the lowest lying states may turn out to be tachyonic. Another 
possibility is that, in the leading order of the large D approximation, the 
instability may not be visible 5 . For example, it is easy to construct a simple 
quantum mechanics problem with D degrees of freedom, where there is a 
metastable state which decays by tunneling 6 . It is usually the case that 
tunneling is suppressed in the leading large D limit, and the metastable 
state becomes stable. One has to go beyond the leading order to see signs 
of instability. It is possible that this is what happens in the model we are 
studying here. 

8. An Additional String Mode 

In this section, we are going to compute a particular correction to the 
leading mean field or large D result 1 . We recall that, in this limit, the 
composite field p can be replaced by its ground state expection value p , but 
to go beyond the leading term, one has to expand in powers of the fluctuations 
around the mean. Of course, an exact computation of the full expansion is 

5 I would like to thank Jeff Grcensitc and Charles Thorn for this suggestion. 
6 I would like to acknowledge a helpful conversation with Eliczer Rabinovici on this 
subject. 

7 See also Appendix B of reference [2]. 
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impractical; however, as in section 6, we look for the dominant contribution 
in the limit of large cutoff. In this case, this is a logarithmically divergent 
term which dominates the rest of the terms, which are finite. To isolate this 
contribution, we split the composite field p into the constant mean value po 
(see eq.(70)), and a fluctuating part x- Along with a power series expansion 
in X) w e treat x as a slowly varying function, and we expand it around a 
fixed point in powers of the coordinates a and r. This latter expansion, 
which is sometimes called the derivative expansion, is very useful in isolating 
divergent terms in the perturbation series and it is widely used in literature. 
The point is that increasing order in this expansion goes with increasing 
number of derivatives on x, and by dimensional arguments, this results in 
greater convergence. The divergent terms therefore appear only in the lowest 
orders of the derivative expansion and they are easy to identify. We shall see 
that the leading divergence is quadratic, but this is already included in the 
calculation done in section 6 with a constant po- The next leading divergence 
is logarithmic, which is the contribution we are going to calculate. The rest 
of the terms in the expansion are finite. The logarithmic term has a special 
physical significance: it provides a kinetic energy term for x i n the action and 
so it promotes x m to a new propagating degree of freedom. This phenomenon 
should be familiar from other two dimensional models, such as the CP(N) 
model [17] or the Gross- Neveu model [18]. In contrast, the finite terms in the 
expansion are non-local and they do not seem to have any special physical 
significance. 

Our starting point is eq.(48) for S m , but with p replaced by p, since we 
are considering fluctuations of p around the mean value po- We define 



is the fluctuating field. Doing the functional integral over q gives the follow- 
ing contribution to the action: 




(70) 



where A is defined by eq.(50) and, 




(71) 




(72) 
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k=0 



Figure 6: First Order Tadpole Graph 

We are going to examine in detail only terms up to second order in x\ it will 
then be easy to figure out the contribution of the higher order terms. We 
therefore expand S' up to second order: 

S' = S' + S' 1 + S' 2 + ..., 
S> = l -DTr\n(dl-A^), 

Si = - % -DA*Tr({dl-A*d^d T xd T ), 

S> 2 = ^DA'Tr^-A^y'drxdridl-A^-'drxdr). (73) 

The zeroth order term S' was called in section 6 and it was already 
calculated there (eq.(50)). The first order term S[ is represented by the graph 
in Fig.6, where the external line in this graph carries zero momentum. Its 
contribution is given by 

S[ — c\ J da J drx((T, t), 

where c\ is a quadratically divergent constant. In addition, there other con- 
tributions from higher order terms, represented by graphs of the form given 
in Fig. 7, which generate the following series in S': 

da drY.c n (x(°,T)) n . (74) 

J 71=1 
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Figure 7: Higher Order Tadpole Graph 



This can be thought of as a potential for x, which, when minimized, deter- 
mines the expectation value of x- However, all this does is to shift the value 
of po, which we have already calculated in the previous section. This ambi- 
guity is due to the intrinsic arbirariness in the split made in eq.(70): Only 
the particular combination of p and x that appears on the right hand side of 
(70) is well defined: Shifting p and the expectation value of x while keeping 
the right hand side of (70) constant will not change anything. We can resolve 
this ambiguity by setting the expectation value of x equal to zero. In that 
case, there is no shift in the value of po, and we can drop the terms given in 
eq.(74). We shall do so in what follows. 

Let us now focus on S' 2 , the term quadratic in x, which is represented by 
the graph in Fig. 8. In momentum space, we can set 

S' 2 = ~DA^ J Ml(%X)x(KX)x(-K,-K), (75) 
where x is the Fourier transform of x, an d 

(4*2 - (k' ) 2 ) 2 



I(ko,h) = J 



d z k- 



OJ 2 



((2fci + k[) 2 - A 2 (2k + k'^) ((2fci - k[f - A\2k Q - k' 

(76) 

So far, this expression is exact, but now, we are going to carry out the 
derivative expansion explained earlier, which coincides with the expansion 
in powers of the momentun k! . The zeroth order term in k' contributes to 
the potential in x, and we have explained above that if one starts with the 
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correct value of the expectation value of p, this term is already taken care 
of. The first order term in k' vanishes, and the second order term has a 
logarithmic divergence. Setting 

I = I + h + h + ■ ■ ■ , 
Ii is zero and I 2 has the form 

h = h fi {Kf + hAK)\ (77) 

and, after Euclidean rotation in k , 

T j2t,^o(3^ 2 ^0^1 + ^l) 

h > 1 - 2J dk (ki + A^kiy (78) 

The integrals are elementary. After the change of variables by 

k\ = r sin(6>), A k = r cos(0), 

the 9 integrals are easily done, leaving a logarithmically divergent r integral: 

in r x dr m f x dr , . 

We have introduced both an infrared cutoff e and an ultraviolet cutoff A to 
regulate the divergent r integral. Putting this back into the equation for S' 2 
(eq.()), and rewriting it in the position space, we have 

S ' 2 ~ ^A J da J dT ( A2 ^ - ^) ' ( 8 °) 




Figure 8: One Loop Contribution to the Propagator 
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where we have defined 




We close this section with a couple of comments on this result: 

a) The contribution of the x mode to the action is exactly the same as the 
contribution of q to S m (eq.(48)), apart from an overall cutoff dependent mul- 
tiplicative constant. This constant can be eliminated by rescaling x (wave 
function renormalization). The important point is that it is positive; other- 
wise, this term would have the wrong (ghostlike) signature. 

b) There is therefore an additional string mode represented by x, with the 
same slope 

a' = A/2 

as the other modes. This changes the dimension of the target space from D 
to D + 1, and therefore it can regarded as a non-leading contribution in the 
large D limit. 

c) If we continued the derivative expansion of S' 2 beyond second order, we 
would get terms with higher derivatives of x with respect to r and a. By 
simple power counting, these terms would be finite and therefore they would 
not be of interest to us. 

d) Let us now consider terms higher then second order in x m the expansion 
of S'. The logarithmically divergent contribution still comes from second 
order in the derivative expansion. Consider the graph of Fig.9, where two 
external lines carry finite momenta and the rest of the external carry zero 
momentum. A simple generalization of the calculation given above in the 
case of the second order term in x shows that this term must be of the form 

S' = JdaJ drh{ X ) {{d TX f - A 2 {d aX f) , (81) 

and h(x) is a function whose power series expansion in x starts with the 
positive constant that appears in eq.(80). Redefining the field x by 

x -> g(x), 

and choosing the function g so that it satisfies 

h(f(x)W(x)) 2 = 1, 

one can get replace h by one. This shows that a simple field redefinition gets 
us back to the second order result given by eq.(80). 



34 



To complete our discussion, we should also consider the contribution com- 
ing from S g in eq.(50). This contribution has a quite different structure 
compared to the one coming from the matter sector. It has a linear cutoff 
dependence, and it depends only on a derivatives of x an d n °t r deriva- 
tives. We remind the reader of the motivation behind integrating over the 
matter fields in S m : The integration produced a kinetic energy term for the 
originally non-propagating field x- Since integrating over the ghost fields 
produces nothing of comparable interest, it is probably best to leave S g as it 
is. 

9. Conclusions 

This article is a direct follow up of reference [2]. The goal is to give a sys- 
tematic treatment which resolves various problems encountered in [2], while 
still staying within the framework of the simple meanfield approximation de- 
veloped in that reference. The crucial improvement over the treatment given 
in [2] is the fixing of an accidental gauge invariance. In addition, there are 
various other technical improvements: Supersymmetry is introduced on the 
world sheet to keep better track of matter-ghost cancellation, and a better 
treatment of the singular determinants is presented. Also, following reference 
[4], we show how to impose partial Lorentz invariance. 

The meanfield method used here and in [2] is sufficiently simple to allow 
the carrying out of a fairly complete treatment. In the leading order of this 
approximation, a condensation of Feynman graphs takes place, which means 
that large order Feynman graphs dominate the perturbation series. As a 
consequence, with a suitable tuning of the coupling constant, a string with 
finite slope is formed. We also show that a new dynamical degree of freedom 
emerges, extending the transverse dimension of the string from D to D + 1. 
Although the string appears to be stable in the leading order, we argue that 
the fundamental instability of the underlying field theory probably shows up 
in the non-leading orders of the approximation. 

The question whether string formation takes place in a given field theory 
is fundamentally important but difficult to answer in practice. The meanfield 
approximation used in this article, free from its earlier shortcomings, is simple 
yet powerful enough to answer this question in the case of the 3 theory. The 
possibility of applying this method to more realistic theories appears more 
appears within reach. Another valuable line of future research would be to 
try to improve over the meanfield approximation. 
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Appendix 

In this appendix, we will discuss an alternative, and more conventional 
way of doing the mean field approximation, and we will show that it is 
completely equivalent to the approach used in this article. Instead of p 
(eq.(39)), we introduce a different set of composite fields: 

AS = J da J dr(/ci(y • q - Ai) + K 2 (h • b - A 2 ) + k 3 (c ■ c - A 3 ) 

+ k 4 (z • r - A 4 ) + K 5 (y • z - A 5 ) + K 6 (b • c - A 6 )) . (82) 

In the large D limit, the fields Ki and A« become classical and they can be 
replaced by their ground state expectation values: 

Ki ->■ (Ki), Xi ->• (Ai). 

The justification for this is well known [16]: The composite fields Aj are each 
sum of D terms, and consequently, they grow like D as D becomes large. On 
the other hand, compared to this, the quantum fluctuations are suppressed 
by a factor of \j\f~D. As a result, in leading order large D these fields 
become classical. In contrast, there is no comparable direct argument for 
why p should become classical in this limit. Therefore, from the perspective 
of large N (large D) physics, the approach sketched in this appendix is better 
motivated than the approach developed in the main body of the article. The 
disadvantage is that, at least initially, many more auxilliary fields have to 
be introduced. We will show presently that the expectation values of these 
extra fields can easily be expressed in terms kq and po introduced earlier. To 
do this, we first write Sf, the fermionic part of the action, in terms of the 
fields introduced in eq.(82): 

Sf = J da J drii^ + -^(4>ii> 2 + foipi) 

+ -(Ai + A 2 + A 3 + p + X 4 - a\ 5 - aX 6 ) ^ + -(A 5 + A 6 )J , 

(83) 



36 



where, to keep things simple, we have set all the /3's equal to one. Comparing 
this with the Sf given by eq.(41) leads to the identification 

= Ai + A 2 + A 3 + p + X 4 - a\ 5 - a\ 6 , (84) 

so that the two expressions for Sf, apart from an additive term, agree. We 
now consider all the A dependent terms in 

AS + S f 

and write down the equations of motion by varying the A's, subject, however, 
to the constraint that (eq.(84)) is held fixed. These equations give 

Kl = K2 = k 3 = — Kb = KQ = a { (85) 

p + a 

Consequently, there is only one independent k, say k%. With the further 
identification 

1 

«i = Po, «i = Po, (86) 

a 

we recover the action of eq.(36) with /3's set equal to one, establishing the 
equivalence of the two approaches. 
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